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z=1 Multifractality of Swift short GRBs? 
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ABSTRACT 

Aims. We analyze and characterize the angular distribution of selected samples of gamma ray bursts (GRBs) from Batse and Swift 
data to confirm that the division in two classes of short- and long-duration GRBs correspond also to the existence of two distinct 
spatial populations. 

Methods. The angular distribution is analyzed by using multifractal analysis and characterized by a multifractal spectrum of dimen- 
sions. Different spectra of dimensions indicate different angular distributions. 

Results. The spectra of dimensions of short and long bursts indicate that the two populations have two different angular distributions. 
Both Swift and BATSE long bursts appear to be homogeneously distributed in the sky with a monofractal distribution. Short GRBs 
follow instead a multifractal distribution for both the two samples. Even if BATSE data may not give a secure interpretation of their 
angular distribution because of the instrumental selection effects that mainly favor the detection of near GRBs, the results from Swift 
short GRBs confirm this behavior, also when are included GRBs corrected by the redshift factor. The distributions traced by short 
GRBs, up to z = 1, depict a universe with a structure similar to that of a disordered porous material with uniformly distributed 
heterogeneous irregular structures, appearing more clustered than what expected. 

Key words. Gamma rays: general, (Cosmology:) large-scale structure of the Universe, Chaos. 



1. Introduction 

Gamma Ray Bursts (GRB) are catastrophic explosions of cos- 
mological origin that illuminate the sky once or twice a day 
with relatively short, intense, flashes of y-rays on the order of 
the MeV and duration that ranges from 10~ 3 to about 10 3 sec- 
onds. Until the launch of Swift, the most widely accepted taxon- 
omy of GRBs is the division between short-hard and long-soft 
bursts, after Dezalay (1992) and Kouveliotou (1993), who found 
in BATSE data a bimodal distribution in the burst duration de- 
fined as the time it takes the 90% of the flux to arrive (Tgo), with 
respect to the local time of the detector. These two distinct GRB 
distributions are separated by a minimum located at Tgo -2 s. 
This sharp division is supposed to be also due to a selection effect 
of the instrumentation onboard the satellite. In fact, new Swift 
observations permit to include in the classification scheme, as 
short bursts, also some distant events having apparently Tyo > 2 
s with respect to Swift's proper time. The determination of the 
redshift permits in fact to correct the Tgo by the relativistic effects 
( |Che et al., 19971 |Ruffini et al„ 2009) . This suggests the need of 
a more solid classification of GRBs on the basis on a broader set 
of criteria, beyond the mere burst duration and the verification of 
the present taxonomy of long/short GR Bs (|Donaghy et al., 2006 , 



Zhang et al., 2007||Bloom et al., 20081 |Belczynski et al., 20 08). 

Long and short GRBs are thought to have been gen- 
erated by different progenitors. Anyway, different progeni- 
tors are thought to produce a black hole accreting material 
from a disc or from a torus with the emission of gamma 
rays (Ghirlanda et al., 2009). The basic model for long-duration 
GRBs, related to the catastrophic release of energy from the col- 
lapse of massive stars (Woolsey, 2001 , Fryer & Kalogera, 2001 1, 
has received a strong support from the observations of their 



X dGehrels et al., 2008| |Nysewander et al„ 2008) , optical and 
radio counterparts, and the association with supernova detec- 
tions, of which the GR B980425/SN1998bw association was 
the first clear exa mple ( [Galama et al., 1998| |Kulkarni, 1998| 
van Paradijs, 1999). Thus, X-ray flashes, associated to long- 
duration GRBs, are probably produced by the highly-relativistic 
jets ejected in core-collapse supernova explosions. The rel- 
ativistic fireball model ([Goodman, 1986 Paczyhski, 1986 



|Paczyhski, 1990|[Castro-Tirado et al., 2001) accounts for a rea 
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sonable description of the observed afterglow spectrum, pow- 
ered by the synchrotron emission of electrons accelerated in a 
relativistic shock with an estimated total energy budget roughly 
the same order of magnitude as that of supernovae Ib/c. The 
discovery of the slowly fading X-ray, from optical and radio 
afterglows of GRBs and the identification of host galaxies at 
cosmological distances, gave further support to describe their 
progenitors as the product of the evolution of short-lived mas- 
sive stars in different cosmological epochs, whose detection 
is limited by the BATSE threshold that sets the observational 
limit at distances slightly larger than z ~ 4 (Wijiers et al., 1997 
Blain & Natarajan, 2000 ). Now the new da ta from Swift passec 
the barrier of z = 8 flSalvaterra et al., 2009) 1. 

Short GRBs are thought to be produced by highly - rela- 
tivistic jets ejected in different processes, as during neutron star 
- neutron star (NS-NS) or black hole - neutron star (BH-NS) 
binary mergers (Narayar i et al., 1992) , whose averaged redshift 
distribution, in the Swift-era, seems to be (z) ~ 1.0 (see e.g. 
Magliocchetti, Ghirlanda & Celotti, 2003, Tanvir et al. 2005, and 
Ghirlanda et al. 2006). Anyway, recent Swift and HETE-2 obser- 
vations put in evidence the presence of clear similarities with the 
afterglows detected in the correspondence of long GRBs, thanks 
to the detection of X-ray and optical afterglows. Other similari- 
ties found between the initial stage of the spectra of short and 
long GRBs together with the presence, in some cases, of X- 
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anism during the first few seconds (Barthelmyet et al., 2005 
|Coward, 2 007). Short GRBs have been observed mostly in el- 
liptical galaxies, but even less frequently also in nearby irreg- 
ular and in star-forming galaxies, confirming as progenitor the 
binary merging scenario (Meszaros, 2006 Narayan et al., 2001 



ray flares, seem to indicate the presence of a common mech- the angular distribution of luminous matter tends to homogene- 
ity. According to Pietr onero, 1987| the distribution of luminous 
matter has multifractal properties with fractal dimensions vary- 
ing from d - 1.23, for the nearby galaxies, up to the value 
d — 2 at very high redshifts. Usually, multifractal distributions 
are present when a structure has different fractal dimensions on 
different parts of the support; in other words, when spatial corre- 
lations are present and change the geometrical shape of the dis- 
tribution at different scales (Falconer, 1990). Those distributions 
cannot be adequately described by a geometrical support with a 
single fractal dimension, but require instead a whole spectrum 
of dimensions. 

To characterize the angular distribution of GRBs in terms of 
multifractals, we use the "method of moments" that estimates 
the fractality of the distribution by calculating the multifrac- 
tal spectrum of generalized dimensions D q in a given range 
( Falconer, 19"90"l|Feder, 19881 Consider N GRBs in a box of size 
L divided into cells of size r. The sample is said to have a non- 
null q - th moment % q on the scale r iff 



Belczynski et al., 2008)) als o through the mechanism of star 
swapping ( |Grindlay, 2006|l. A similar b ehavior is expected 
for NS-BH binaries (Paczyhski, 1991 1. The angular and 
spatial distributions of BATSE GRBs look quite isotropic 
(Briggs et al., 1996), with only very small anisotropics observed 
in the angular distribution (Meszaros et al., 1999| ); short- and 
long- GRBs in BATSE catalog actually show two different an- 
gular distributions ( |Vavrek et aT72 008 Bal azs etaL72 009), but 
the connection to the hypothesis of the instrumental selection 
effect is not still completely clear. 

In this paper we characterize the angular distributions of the 
two classes of GRBs with multifractal analysis from a selected 
sample of BATSE and Swift observations, taking also in con- 
sideration the redshift relativistic effects and compare the results 
obtained from the two catalogues. In Chapter 2 are described the 
mathematical basis of the method. In Chapter 3 are estimated 
the fractal/multifractal dimensions by determining the distribu- 
tion moments of the multifractal spectrum from the second up to 
the tenth order and then draw the conclusions. 



2. Fractal / Multifractal analysis of the angular 
distribution 

Being GRBs associated to galaxies, they should trace the angu- 
lar distribution of their host galaxies at distances slightly larger 
than those estimated with classical Supernovae, i.e. distances 
the order the Gpc ( [Paczynski, 1986| |Usov & Chibisov, 1973) , 
across which inhomogeneity in the distribution of luminous mat- 
ter should be averaged out from the Mpc scale. 

Mandelbrot and Peebles (1980) first used fractal geome- 
try to describe the angular and spatial distribution of galax- 
ies. In their work they used Levy-Rayleigh random paths, i.e. 
infinite-variance, stable, generalized random walks in which the 
step lengths are described by a tailed probability distribution. 
Galaxies are placed at the steps of a random walk and each 
galaxy is randomly connected with another in the close proxim- 
ity, like mimicking the random motion of a fly in the air passing 
through each given point. The power-law constraint of the mo- 
tion gives the fractal exponent of the random path and rules the 
distribution of the jump lengths and the direction of each jump 
is taken isotropically at random (Martinez, 2002). 

Galaxy surveys show that the distribution of the luminous 
matter in the universe is more complicated than that of a single 
fractal. The distribution of galaxies has instead multifractal prop- 
erties (jPjetronero, 1987 , |Celerier & Th ieberger, 200T) , with a 
tendency to cluster following the well known peculiar structures 
made by patterns of voids and filaments (Kurokawa et al., 1999 ). 
The CfA survey of nearby galaxies, for example, presents scal- 
ing properties with a correlation dimension in the interval D2 — 
1.3-1 .4. The two-point correlation function and the power spec- 
trum analysis indicate that the distribution of galaxies at very 
large scales becomes homogeneous and isotropic like the X-ray 
background emitted by active galactic nuclei (AGNs), in agree- 
ment with the Cosmological Principle (Peebl es^ 1993| ). 

Multifractal scaling analysis revealed the presence of ag- 
gregated structures also in some samples of galaxies with dis- 
tances larger than r ~ 20 Mpc. At scales larger than 30 Mpc, 
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where n, is the number of objects present in the i-th cell. 
Higher moments reflect the emergence of structures present in 
the denser regions. From the scaling relation conducted on a two- 
dimensional section we obtain (Kurokawa et al. 2001), 
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where the coefficients D q belong to the spectrum of generalized 
dimensions. Do is the capacity dimension, D\ is the information 
dimension, and Di is the correlation dimension. 

To estimate the multifractality of the distribution one has 
to determine each generalized dimension D - q as a function 
of the q-th moment. One gives evidence to the dependence of 
each of the fractal dimensions from the moments by drawing 
the q - D q plot, whic h is based on the Lipshitz-Holder exponent 
dBenassi et al., 1998|[Goltz, \991) . In the q-D q plot one defines 
monofractals those simple fractal structures described by only 
one dimension, Df, and in the plot all the D 9 's result equal to 
Df. In the general case of a multifractal, instead, D q is usu- 
ally decreasing for higher and higher values of the moments q: if 
q > q', then D q < D q > until converging to the asymptotic value 
Doo for a distribution of infinite objects. 

To characterise with a good precision the multifractal distri- 
bution of a sample of object in space one usually does not need 
to calculate the generalized dimensions up to the limit dictated 
by the number of points in the space. For our purposes a good 
estimation was obtained by taking as upper limit q — 10. 

3. Results and discussion 

From Swift recent observations we selected a sample of 444 
GRBs without ambiguity of classification, 52 of which are clas- 
sified as short- GRBs after the correction of T90 from the rela- 
tivistic effects (data updated at 2009-09-15). From BATSE 4 cat- 
alogue, instead, the total number of GRBs used in our analysis 
is 1843, 1447 of which are classified as long-GRBs. 

To verify whether the classification in the two subgroups ac- 
tually correspond to two distinct populations in both the sets of 
data collected by Batse and Swift, we also performed a series 
of tests of the angular distributions in different sub-samples, by 
randomly choosing from our data a sets of either short or long 
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GRBs and mixing them together, to see whether GRBs might be 
discriminated without choosing "a priori" the two classes fol- 
lowing the burst duration time distribution. The two classes of 
long/short GRBs emerged anyway from the fractal analysis of 
their angular distribution in the sky. More precisely, by progres- 
sively mixing the population of short GRBs with a growing sam- 
ple of randomly chosen long GRBs, the multifractal dimension 
converges to that of a homogeneously distributed monofractal 
with dimension d — 2. This simple test corroborates the ac- 
tual sub-division of GRBs in the two short/long populations for 
both the Batse and Swift data and also when the two datasets are 
mixed together. 

We now present and discuss the results obtained for each 
different class of GRBs: 

Long bursts: Fig. (Q]) clearly shows the uniform angular dis- 
tribution of both the samples of Swift and BATSE 4 selected 
data. The sample of data from Swift shows a homogeneous dis- 
tribution. In this case, the moments D q tend to decrease slowly 



Swift Long GRBs q - Dq plot 
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Fig. 1. Upper panel: q — D q plot of Swift long GRBs. The long 
bursts follow an angular homogeneous distribution. The mo- 
ments D q tend to decrease slowly from D\ ~ 2 to D q ~ 1.8 with 
a rise of the fractal dimension after D q = 8. The spectrum of 
dimensions is that of a two dimensional brownian motion with 
fractal dimension 2. Lower panel: also BATSE 4 long bursts 
present a very similar homogeneous structure, with a mean frac- 
tal dimension D 1.98. In this case, a deviation in the distribu- 
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creasing after the moment q = 8. In the plot are also present 
some sporadic and discontinuous jumps down to D q ~ 1.7 that 
deviate from the main smoothness, a behavior that might be due 
to numerical errors in the determination of certain D q 's . In any tion due to the moments 4 < D q < 6, is present, 
case, this does not give rise to misunderstandings to the interpre- 
tation of the global behavior observed in the plot. What is seen 
is a typical example of a stochastic homogeneous distribution, 
similar to that of a fractional brownian motion (FBM) having di- 
mension D — 2, as expected from distant sources that homoge- 
neously distribute according to the cosmological principle. The 
errorbars of the D 9 s' are calculated from the error propagation of 
the statistical uncertainty of the position of each GRB and from 
the instrumental errors. In the case of Swift data, errorbars are 
too small to be visible in the graph. 

By analyzing BATSE data, we find a different behavior in 
the interval 4 < D q < 6, where is seen a small bump. Also in 
this case, within the experimental errors, the angular distribution 
is characterized by an almost constant fractal dimension fluc- 
tuating within the interval D q = [1.8 - 2] that approximately 
describes the angular distribution of a uniform structure. Having 
found a similar result in both the datasets, being the redshift scale 
around the Gyr (as confirmed by the redshift values obtained for 
most of the bursts), long GRBs distribute according to the the 
Cosmological Principle, in which the Universe is completely ho- 
mogeneous and isotropic at large scales. 

Short bursts actually represent a more articulated and com- 
plicated set of phenomena with still unknown aspects. Fig. [2] 
shows that the whole sample of 53 Swift short bursts (together 
with the apparently-long bursts then corrected by z and those 
with Tcjo already below 2 s) is characterised by a spectrum of di- 
mensions that decreases for increasing g's, which is the signature 
of a multifractal distribution. 

Short GRBs look angularly distributed as they were driven 
by a multifractional brownian motion (mBM) of mean index 
H = 0.68, a distribution that starts with a clustered shape then 
evolving towards homogeneity at larger scales. The mBM is an 
extension of the fractional Brownian motion in the sense that 
the path regularity can vary with time as is seen in the evolu- 
tion of anomalous diffusion processes. The Universe traced by 
short short GRBs exhibits a multifractal structure, correlated at 
lower scales and evolving to homogeneity at larger distances. 
This structure appears to be similar to a disordered porous mate- 
rial, that presents heterogeneous structures, or even irregular in a 
uniform sense. This type of foam contains multiple, nested natu- 
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Fig. 2.q-D q plot of Swift short GRBs. Upper panel: spectrum 
of fractal dimensions of the angular distribution of the whole 
sample. This behavior with a maximum around D q = 3 clearly 
indicates a multifractal homogeneous distribution very close to 
that of a multifractional brownian motion. Lower left panel: di- 
mension moments of the z-corrected GRBs with relatively high 
redshift. The structure is clearly multifractal. Lower panel: also 
those GRBs with already < 2 sec. show a multifractal dis- 
tribution, but with an asymptotic convergence to D q = 1. In this 
case the structure traced by the GRBs appears more clustered. 



ral length scales or continuously evolving scales, while moving 



to higher redshifts ( |Lim and Muniandy, 2002 1. 

The further division in two sub-classes with respect to the 
rcjo-corrected by the redshift and those that were already below 
the 2 seconds, show that the two sub-samples follow two mutu- 
ally different multifractal distributions and the sub-class of the 
z-corrected short GRBs presents a higher fractal dimension qo 
which is close to that of a homogeneous distribution, confirming 
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Batse short bursts: q - Dq plot 
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Batse Short (Class III) q - Dq plot 
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Fig. 3. Upper panel: q - D q plot of BATSE 4 short GRBs (Class 
I + Class II) present a fractal structure with dimension fluctuat- 
ing around the value > 1.8. Lower left panel: q - D q plots of 
BATSE 4 short bursts Class II and Class III alone clearly follow 
a multifractal distribution. 



the behavior expected from a mBM (see the two lower panels in 
Fig.E). 

The whole sample of BATSE short GRBs shows an articu- 
lated spectrum of the general dimension D q , with a smoother be- 
havior with respect to the whole sample of BATSE long GRBs, 
as reported in the upper panel of Fig. [3] The capacity dimen- 
sion Dq and the correlation dimension Do present similar values 
close to D q =i 1.9 that are smaller than those of long GRBs. 
This demonstrates that short GRBs have a different multifrac- 
tal and more clustered distribution than the class of long ones. 
This difference is better put in evidence if BATSE GRBs are di- 



vided in the two subclasses indicated bylMukherjee et al., 1998 



Class II with short/faint/hard bursts and Class III with intermedi- 
ate/intermediate/soft bursts (we remind that Class I correspond 
long GRBs), as drawn in the lower panels of Fig [3] Here, the 
generalized dimension D q of both II and III classes alone indi- 
cates that the angular distribution is clearly multifractal, similar 
to that seen in the galaxy distribution at low redshifts. 

Concluding, even if this effect could be ascribed to a selec- 
tion effect in BATSE 4 data, anyway Swift results confirm the 
existence of two different angular distributions associated to the 
two classes of short and long GRBs. In both Swift and BATSE 
data sets long bursts are homogeneously angularly distributed. 
Short GRBs, instead, trace a distribution which appears differ- 
ent to that expected from the clustering of luminous matter in 
the Universe around the redshift value z ~ 1. We in fact still 
clearly see a multifractal distribution in a structure that should 
be already homogeneously-distributed. By taking in account the 
fact that the candidate progenitors of short GRBs are thought 
to migrate away from their initial positions, we should also ex- 
pect an additional blurring effect on the detected structure that 
should smear out the distribution towards the total homogeneity, 
which has not been observed. This suggests that either GRBs are 
not good tracers for the matter distribution or that the Universe 
traced by GRBs might appear more clustered at redshifts z < 1 
than expected. 
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